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Abstract 


We study the dependence of the exclusive decay mode B^DtV in type 
II two Higgs doublet models and show that this mode may be used to put 
stringent bounds on tan/3/m//. There are currently rather large theoretical 
uncertainties in the distribution, but these may be significantly reduced 
by future measurements of the analogous distribution for B^D{e, p)v. We 
estimate that this reduction in the theoretical uncertainties would eventually 
(i.e., with sufficient data) allow one to push the upper bound on tan/3/m// 
down to about 0.06 GeV“^. This would represent an improvement on the 
current bound by about a factor of 7. We then apply the method of optimized 
observables which allows us to estimate the reach of an experiment with a 
given number of events. We thus find that an experiment with, for example, 
10^ events could set a 2a upper bound on tan fl/mn of 0.07 GeV“^ or could 
differentiate at the 4.6(T level between a 2HDM with tan/3/m// = O.I GeV“^ 
and the SM. 
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I. INTRODUCTION 


There has recently been considerable interest in constraining the parameter space of 
type II two Higgs doublet models. The main reason for this interest, of course, is that the 
Higgs sectors of minimal supersymmetric extensions of the standard model are generically 
of this type |I|]. The charged Higgs sectors of two Higgs doublet models (2HDM’s) may 
be characterized by the ratio of the two Higgs’ vacuum expectation values, tan/3, and the 
mass of the charged Higgs, niH- In this work we will investigate how the exclusive decay 
channel B DtV may be used to place tight constraints on the ratio tan/3/mj:^. This 
channel is expected to have a branching ratio on the order of half a percent |Q, so that one 
would expect on the order of 10® such decays at the B factories which are currently under 
construction. 

There already exist several constraints on tan (3 and mn- The most direct lower bound on 
the charged Higgs mass comes from the non-observation of charged Higgs pairs in Z decays 
and gives tuh > 44 GeV [Q. Another limit comes from top decays, which yield the bound 
rriH > 147 GeV for large tan/3 0]. Finally, for pure type H 2HDM’s one finds mn > 300 
GeV, coming from the virtual Higgs contributions to 6 —> sy [Q . This latter limit disappears 
in the context of supersymmetry since the Higgs contributions to 6 —sy can be cancelled 
by other contributions [^. There are no experimental upper bounds on the mass of the 
charged Higgs, but one generally expects to have mn < 1 TeV in order that perturbation 
theory remain valid [j^. A lower limit may be placed on tan/3 by considering the branching 
ratio for Z —> bb. The resulting bound of tan/3 > 0.7, obtained in Ref. [||, coincides with 
the range generally favoured by theorists in order that renormalization group evolution drive 
electroweak symmetry breaking [W|. For large tan/3 the most stringent constraints on tan/3 
and rriH are actually on their ratio, tan (3 /rriH- The current limits come from the measured 
branching ratio for the inclusive decay B —> XtV, giving tan(3/mH < 0.46 GeV“^ 
from the upper limit on the branching ratio for B 


11 , and 


rz/, giving tan/3/mH < 0.38 GeV ^ [jT2 


While both of these limits are quoted as being at the 90% confidence level, the latter may 
be somewhat less constraining due to the uncertainties in 14 b and fs- 

Our main goal in this paper is to investigate the sensitivity of the exclusive decay B 
DtV to the ratio tan j3 /rriH- We will concentrate on the distribution for this decay 
and discuss how the theoretical uncertainties in this distribution may be minimized. We 
also apply the optimized weighting procedure ||l3|,p!^ to this distribution in order to derive 


quantitative estimates for the sensitivities of experiments with given numbers of events. 
This procedure can be shown to give the smallest statistical uncertainty when analyzing 
the data in a given experiment. The present work complements the previous theoretical 
studies of the inclusive 


T5|-^ and exclusive |E3 


of the purely leptonic decays B ^ tv [0 and B^ ^ tv 


semi-tauonic B decays, as well as those 
These decays are attractive 


because the Higgs contribution occurs at tree-level and cannot be cancelled by, for example, 
supersymmetric loop effects. Thus, the results of our analysis should be applicable to any 
type H 2HDM This situation may be contrasted with that in 6 —>■ sy [Q. 

A feature which is common to all of the tauonic and semi-tauonic B decays is that the 
Higgs contribution to the amplitude interferes destructively with that due to the standard 
model (SM). As a result, the corresponding integrated partial widths, plotted as functions 
of tan 13/mu, tend to have minima around tan 13/mu ~ 0.2 - 0.3 GeV“^. Most studies 
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to date have concentrated on the region to the right of the minimum, where the Higgs 
contribution to the width begins to dominate over the SM contribution. Indeed, the present 
experimental limits - derived using integrated partial widths - correspond to this region. 
In order to use semi-tauonic decays to probe values of tan/3/m// near and/or below the 
minimum, it will be extremely useful to have detailed theoretical predictions for quantities 
beyond simply the integrated partial widths. This is because the plots of the widths as 
functions of tan/3/m// are generically relatively flat up to tan/3/m// ~ 0.4 GeV“^. Several 
authors have suggested using the energy distribution or longitudinal polarization of the r 
in this regard, but this may be difficult experimentally since two neutrinos are always lost. 
An alternative approach, which we will study in detail, is to use the distribution. A 
possible drawback of this approach is that the distribution is very sensitive to theoretical 
uncertainties in the shapes of the hadronic form factors. As we shall see, however, the 
situation in the exclusive channel B DrV appears to be quite encouraging. The reason 
for this is that once the distribution for B D{e, ii)V has been measured, that for B —> DtV 
can be predicted with relatively small theoretical uncertainties. The resulting distribution is 
quite sensitive to tan (d/mn, even for relatively small values of this ratio. Furthermore, the 
distribution has a qualitatively different shape for values of tan/3/m// above and below 
the “critical value”, tan/3/m// ~ 0.3 GeV“h 

We have chosen to focus on the decay channel B —> DrV instead of on H D*tV, even 
though the latter channel will likely have a somewhat larger branching ratio and may also 
be more accessible experimentally (in analogy with the decays to the lighter leptons p7| ) 
than the former. Our main motivation for considering B DtV rather than B — D*tV 
is simply that the Higgs contribution has a much larger effect in the former case. This 
feature has already been noted in Ref. and is in part due to an enhancement by a 
factor (m^ + m//)/(mB — m/j) ~ 2 in the effective interaction. As noted in Ref. ||24|| , 
this enhancement effect means that the exclusive D channel is also more sensitive than the 
inclusive channel, since the less-sensitive D* mode tends to dilute the inclusive measurement. 

The plan of the remainder of this paper is as follows. We begin in Sec. || by deriving the 
g^ distribution for B DtV in terms of the dimensionless variable t = q'^/m%. In Sec. |T| 
we estimate the theoretical uncertainties in this distribution and in the integrated width 
once the distribution for B Z3(e,/x)r' has been measured. Barring any further input, 
these uncertainties would eventually limit the reach of such an experiment. In Sec. ITV 


we 


apply the optimized weighting procedure to the g^ distribution and in Sec. 0 we present our 
conclusions. 


II. CALCULATION OF THE DIFFERENTIAL DISTRIBUTION 

The two diagrams which contribute to the decay B^DtV in a type H 2HDM are shown 
in Fig. |I|. The amplitude corresponding to the SM hF-exchange diagram (Fig. |I|(a)) is given 
by 


ATsM = -2V2GFVcb{D{p')\cLYbL\B{p))TL{Pr)'lti^L{Piy), (1) 

where 'i/L=|(l—7^)'0- The matrix element of the axial vector current in the above expression 
is identically zero since one cannot form an axial vector using only p and p'. The vector 
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current matrix element may be expressed in terms of two form factors, Fq and Fi, which are 
dehned as follows: 


{D{p’)\cj>^b\B{p)) = F,{t) 


{p + pT - 


m%-m% 


+ Fo{t) 




( 2 ) 


with q = p — p' and t = /m\. The form factors Fq and Fi are normalized such that 

Fo(0)=Fi(0). There is thus no singularity at g^=0. 

The parametrization of the form factors given in Eq. is particularly well-suited for 
our purposes since Fo(f) and Fi[t) may be associated with the spin-0 and spin-1 components 
of the exchange particles, respectively [^. The contribution to the total amplitude coming 
from the (spin-0) charged Higgs diagram (Fig. |^(b)) may then be included by the following 
replacement in the SM expression for the amplitude: 


F,{t)^Fo{t){l + 5H{t)) 


(3) 


The function Suit) is given by 


bnit) 


/ tan/? y .2 Fs{t) 

\mH ) {ms - ttid) V mb ^) Fo(f) ’ 


where the scalar form factor Fs{t) is dehned by 


{D{p')\ch\B{p)) = {mB + mD)Fs{t). 


(4) 

(5) 


It is now straightforward to work out the expression for the differential partial width in 
terms of these form factors. Let us hrst dehne the following dimensionless quantities: 


td 


ml 
2 ’ 
ml 



The expression for the width is then 

dT{B ^ DtV) _ Gl\V,b\^m% 
dt 1287r3 


( 6 ) 

(7) 


where the dimensionless Dalitz density, p{t), may be decomposed into spin-0 and spin-1 
contributions as follows, 

p{t) = + Po{t) + Pi{t), ( 8 ) 

with 

Poit) = Y [Foit)f (9) 

Piit) = ‘^[Fi{t)f (l-j^ + X^l,rD,t), ( 10 ) 

and 
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A(ci, b, — Qi b ~\~ c — 2(ci6 -I- cic -|- bc^. 


( 11 ) 


The above expression for the differential width is a sum of two semi-positive definite 
terms, corresponding separately to the spin-0 and spin-1 contributions. The spin-0 contri¬ 
bution disappears in the limit Tt-—> 0, so that in this limit we recover the familiar expression 
for the semileptonic decay to an electron or muon. This observation is actually very impor¬ 
tant, since the distribution dV(^e^^)/dt is expected to be measured very precisely at the B 
factories which are currently under construction. Such a measurement would yield valuable 
information regarding the distribution since 


dT{T) 

dt 


dt \ t 



+ 


3r^ (1 - rpf 
2t X{l,rD,t) 




where r^- < t < (^1 



and 


( 12 ) 


CW - (13) 

Thus, the measurement of the differential distribution for the decays B D{e,ij)T7 may be 
used to predict the SM distribution for the decay into the r, up to the function As 

we shall see in the next section, C,{t) may be calculated within the context of Heavy Quark 
Effective Theory with a relatively small uncertainty. 

In the remainder of this work we examine the t distribution in Eq. (|^) in detail and 
evaluate its sensitivity to a Higgs signal. 


III. THEORETICAL UNCERTAINTIES IN THE RATE AND IN THE 
DIFFERENTIAL DISTRIBUTION 

Our goal in this section is to estimate the theoretical uncertainties in the distribution 
dr( t) /dt and in the integrated width if the analogous distribution dT (e^^) /dt has been mea¬ 
sured very accurately. The first step is to evaluate the form factors Fq, Fi and Fs in the 
heavy quark symmetry limit using the results of Heavy Quark Effective Theory (HQET). 
Since in this limit all of the form factors have a known dependence on a single universal 
function - the Isgur-Wise function [^] - one is potentially in very good shape for trying 
to disentangle the Higgs contribution from the SM contribution in the t distribution. The 
symmetry-breaking corrections to this picture introduce theoretical uncertainties into the 
calculation of the function ({t) = To(t)/Ti(t). It is this function which determines the shape 
of the distribution dVi^T-^/dt and whose uncertainties we shall need to estimate. 

It is natural in the context of HQET to define quantities in terms of the meson velocities 
instead of in terms of their momenta. The matrix element of the hadronic vector current is 
then usually written as 


{D{v')\c'y^b\B{v)) = ^/mBmo [h+(tc)(n -f v'Y + h_{w){v — v')^] 


(14) 


where v^=p^/rriB-, v'^=p '^/ttib and 
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, I + Td -t 
W = V ■ V = -- 


(15) 


Comparing the expressions in Eqs. (EP and (|^), we find the (exact) correspondence 




2 r 


1/4 

D 


^ 1 + ^ 


h+{w) - 


. 1 - 


h_{w) 


Flit) = —ui 1(1 + VrD)h+M - (1 - ^/^)h.[w)] 

2 I'D 


(16) 

(17) 


The formalism of HQET gives a self-consistent way to express hadronic form factors 
in an expansion in powers of A/tuq, where rriQ represents the masses of the heavy quarks 
involved in the transition and where A represents a dimensionful quantity which is generically 
of order Aqcd- For the meson form factors the first-order corrections are proportional to 
A/mg, where A is defined as the difference between the meson and quark masses^: 


A = Mm - mg. 


(18) 


To leading order in the 1/mg expansion all of the form factors may be expressed in terms of 
the universal Isgur-Wise function, .^(tn), which satisfies the normalization condition ,^(1)=1. 
In the heavy quark symmetry limit, and ignoring short-distance QCD corrections, one finds 
that h+{w)^^{w) and h_{w)^0, so that 


FM 

Fdt) 


HQS 


1 ft - {1 + y/rBf 


2 r 


1/4 

D 


1 + y/n/ 


HQS 1 + y/n/ 


2 r 


1/4 

D 






Similar considerations for the scalar matrix element yield 


(19) 

( 20 ) 


Fs{t) 


— F„(«) 


( 21 ) 


The above expressions receive corrections due to the finite masses of the heavy quarks. The 
1/mg corrections to h± have been considered in detail in Ref. [^, while the corrections to 


the scalar matrix element do not appear to have been calculated. For this reason we will, 
for the purpose of estimating the theoretical errors, take the relation in Eq. ( 0 ) to be exact. 
This will not lead to significant errors in attempting to bound small values of ta.n j3/mi{. 
Under this assumption, 6nit) takes the simple form 


dnit) ~ 


/tan/3\^ m^m'^^t 

\ mn J {ms - moY 


( 22 ) 


^Note that the pseudoscalar and vector mesons corresponding to a given heavy quark are degen¬ 
erate in mass at zeroth order in the heavy quark expansion. Mm is thus not the mass of any 
particular physical meson j^]. 
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where we have also dropped the term proportional to rric, since its contribution to the 
amplitude is typically very small for the range of tan (3 which we will be considering. 

The form factors h± receive both short- and long-distance corrections. The short-distance 
corrections are embodied in the Wilson coefficients and may be calculated reliably using 
perturbation theory and renormalization group evolution. They typically give corrections 
to the tree-level results which are on the order of 10% [^]. The long-distance corrections 


are intrinsically non-perturbative and give rise to new sub-leading universal functions. At 
order l/mg there are four such functions [^. The corrections to the results obtained in 


the heavy quark symmetry limit may be taken into account by writing 




(23) 


where the renormalized Isgur-Wise function is dehned such that it still satishes 'Cren(l)= 1- 
The explicit expressions for N± to order l/mq may be found in Ref. [^. Their values at 
zero recoil (r(;=l) are of particular interest since one may use this information to extract 
Vcb from the differential distribution for B —>■ D{e,ii)v. It is known that is protected 
by Luke’s theorem [^] and thus does not receive any l/mg corrections at zero recoil. Re¬ 
does receive l/mg corrections at zero recoil, but the resulting contributions to the decay 
rate are parametrically suppressed by the factor [(m^ — mr,)/{mB + ~ 0.23 ||3T|-|33|. 

The uncertainties in the decay rate at zero recoil due to N_ are thus on the order of a few 
percent, which is about the same size as the expected I/uiq corrections. 

We are now in a position to calculate the function C(f)) fhe ratio of Fq and Fi, which 
appears in the expression for the differential distribution given in Eq. (|I2[) . Since |A^_|-C|A^+|, 
as will be clear a posteriori, it is an excellent approximation to expand ({t) in powers of 
N_/N_^.. To hrst order this gives 


C{t) ^ Coo(t) + 


N_{w) 

{l-rD){l + V^yN+{wy 


where 


Coo(t) 


(1 + 


(24) 


(25) 


is the result obtained in the limit of heavy quark symmetry. 

The theoretical error associated with ((t) is actually remarkably small. The main reason 
for this is that there is an approximate accidental cancellation in the expression for N_ 
which tends to make it a very small number, to leading order in l/rriq |^. The resulting 
uncertainties due to our ignorance of the forms of the sub-leading universal functions tend 
then also to be small (on the order of several percent). By way of contrast, the uncertainties 
in are rather large (on the order of 20%), but these have a negligible effect on the 
ratio N_/N^. This means that while the current predictions for Fo{t) and Fi{t) have large 
theoretical uncertainties, the prediction for the ratio of Fo{t) and Fi(t) has a relatively small 
theoretical uncertainty. Thus, once Fi{t) has been determined experimentally, F^t) is also 
known quite well. The ratio N_/N_^ may be estimated by using the forms predicted for the 
sub-leading functions in specihc model calculations. These functions have been calculated 
using QCD sum rules in Refs. We have studied N_/N^ numerically by allowing 
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the sub-leading functions to vary over the regions suggested in the plotsQ in Ref. and 
by allowing A to vary in the range 


0.4 GeV < A < 0.6 GeV. 


(26) 


We have, for consistency, taken the heavy quark masses to be dehned in terms of A through 
Eq. ([T^), setting rrif, = 4.8 GeV and rric = 1.45 GeV when A = 0.5 GeV. The resulting range 
for V_/V+ is then given by —0.06 < N_{w)/N+{w) < 0.0, for 1.0 < tc < 1.6. In order to 
account for the unknown l/rriq corrections, we conservatively double the size of this region. 
We thus estimate the range of theoretical uncertainty in to be 


-0.09 < 


N_{w) 

N+{w) 


< 0.03. 


(27) 


It is now straightforward to use Eq. (^) to determine the differential width for B DtV 
once that for B ^ D{e, ij)V is known. The resulting curve will have a theoretical uncertainty 
determined by the uncertainty in the ratio V_/V+. We illustrate this procedure in Fig. 2 by 
plotting (1 -|- SHit))'^po{t) (shaded bands) and pi{t) (solid line) for the SM and for 2HDM’s 
with tan(3/mH = 0.06, 0.25, and 0.35 GeV“^. The spin-1 contribution to the width, pi{t), is 
independent of tan (3/mH and is assumed to have been determined experimentally from the 
decays to the lighter leptons. The spin-0 contribution, (1 -|- 6H{t))‘^po{t), is determined from 
the spin-1 contribution by using the relation Fo(t) = ((t)Fi(t), with ({t) as given in Eq. (p4|). 
For the purposes of this plot we have used the simple heavy quark symmetry relation for 
Fi(t) given in Eq. (^OD, taking^ ^(tc) = 1 — 0.75 x (tc — 1). Note that the spin-0 contribution 
is qualitatively quite different for tan (3/mH = 0.25 GeV“^ and tan (3/mH = 0.35 GeV“h 
These values fall on either side of the “critical value”, tan (3/mH ~ 0.3 GeV“^, for which the 
integrated width is at a minimum. From this plot we estimate that the current theoretical 
uncertainty would allow one to use the differential distribution in R —> DtV to rule out a 
2IIDM with tan (3/mH > 0.06 GeV“^. 

It is also useful to examine the behaviour of the integrated width as a function of 
tan (3/mH- This behaviour is illustrated in Fig. ^ in terms of the dimensionless quan¬ 
tity p ^ J* p(f)df, which is normalized by the factor Vcfo|^r?7.^/1287r^ (see Eq. (0))- 
As in Fig. H, the shaded band corresponds to the theoretical uncertainty in the ratio 
(({t) = Fo(t)/Fi(t) and does not include the uncertainty in the form factor Fs(t). For 
the sake of illustration we have again used the simple form for Fi{t) given in Eq. (pOD, 
taking ^{w) = 1 — 0.75 x {w — 1). The origin in this plot (i.e., the point tan (3/mH = 0) 
corresponds to the SM. For tan (3/mH > 0, the Higgs contribution begins to intefere with 
the SM contribution, leading at hrst to a reduction in the width and then, for large values of 
tan (3/mH, to an enhancement. For tan (3/mH ^ 0.45 GeV“^, the Higgs contribution com¬ 
pletely dominates the width. It is from this region that the current inclusive semi-tauonic 


^We have allowed the sub-leading function xi®'^(rc) to vary over the larger shaded region shown 
in Fig. 5.5 in Ref. |27|. 


^This form for (,{w) is consistent with the current experimental situation [Q. 








bound on the ratio comes. One may, in fact, compare our plot with the analogous plot 
in the inclusive case (see, for example. Fig. 1 in Ref. |^^). Such a comparison shows that 
while the two plots are qualitatively similar, the curve in the present case has a more pro¬ 
nounced dip near ta.n jS/mn ~ 0.3 GeV“^, dropping to about 50% of the SM value at this 
point. In the inclusive case the branching ratio at the minimum drops to about 80 — 90% 
of the SM value. This feature illustrates a trend which we have already mentioned above: 
since the B —>■ D*tV channel is not very sensitive to the Higgs and since it has a relatively 
large branching ratio, it tends to “dilute” the inclusive mode, making it less sensitive to the 
Higgs contribution. Fig. ^ also illustrates why it is useful to have additional information 
besides the integrated width. If, for example, the measured width is near or below the SM 
value, there will alway be a two-fold degeneracy in the corresponding value of tan/^/m^- 
The differential distribution may be used to differentiate between the two values, however, 
since this distribution is qualitatively different for values of tan 13/mn above and below the 
critical value. The distribution could also be extremely useful in ruling out small values 
of tan/9/m//, since the integrated width is quite flat in this region. 

The shaded bands shown in Figs. and ^ correspond to our estimates of the current 
theoretical uncertainties in Fo(^) and F\(t) and do not take into account the uncertainties 
associated with Fs{t). These uncertainties could in principle be on the order of 10 — 20%, 
although they could also be small, as was the case for N_. It is beyond the scope of this 
paper to provide a more quantitative estimate for the uncertainty in Fsit). Let us note, 
however, that these uncertainties have a negligible effect for small values of tan/3/m// and 
thus do not affect one’s ability to rule out a 2HDM with a small value of tan/3/m//. Should 
one observe evidence for a non-zero value of the ratio tan/3/m//, one would clearly want to 
calculate Fs{t) more carefully in order to precisely determine this ratio. 


IV. USING THE OPTIMIZED WEIGHTING PROCEDURE 

We have so far considered the theoretical uncertainties which arise in the calculation of 
the differential distribution for B DtV. These uncertainties determine - in the limit of 
inhnite experimental statistics - our ability to distinguish the standard model from a two 
Higgs doublet model with a given value of tan/3/m//. Let us now turn the situation around 
and ask the following question: Suppose all of the form factors could be determined precisely 
by some means[|. How well could one then differentiate between the SM and a 2HDM with 
a hnite number of events? This question may be answered by using the optimized weighting 
procedure which has recently been discussed in Refs. P!3|JT^ . This procedure provides the 
most efficient way (with regard to statistical uncertainties) to analyze the experimental data 
in order to differentiate between the two model^. 


^ This is not an unreasonable assumption, since the form factors can in principle be determined 
on the lattice |^]. Even without lattice results, the sub-leading universal functions may eventually 
be constrained by precision measurements in some of the other B decay channels. 

® Note that, while this technique provides the best way to minimize statistical uncertainties, in 
an actual experimental setting one will have to consider the effect of systematic errors as well. 
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In this section we will briefly review the optimal observables method and then apply it 
to the problem at hand, which is to distinguish between the SM and a 2HDM. Let us first 
review the method. Suppose one has a distribution of the form 




(28) 


where (j) represents some collection of kinematical variables, the functions fi{4>) are known 
functions of those variables and the constants Cj parametrize the different models which 
one wishes to investigate. It may be demonstrated that the optimal observables method 
provides the most efficient way in which to extract the coefficients c*. This technique was 
first discussed in Ref. |^] for the case in which Ci = 1 and C 2 = A, with A being some small 

since for large values of 


14 


number. We will use the generalized version presented in Ref. 
tan 13 /thh the q need not be small. 

The main goal of the optimal observables approach is to And the optimal set of functions 
Wi{(()) such that 


Ci = y Wi{(j))0{(j))d(j). 

As shown in Ref. |^, this set is given by 




0{<P) 


where 


= M^\ 


M,, = 


_ f /»(0)/i(0) 

0 ( 0 ) 


dcf). 


The coefficients c* are then given by 

Ci = (y . 


(29) 

(30) 

(31) 

(32) 

(33) 


The experimental task reduces to measuring the elements of the matrix M. The statistical 
error associated with this procedure is embodied in the function, defined by 

= — I](ci - c°)(cj - c5)My, (34) 

ij 

where the represent the measured values of the coefficients, N is the total number of 
events and ar = f 0{(j))d(j). 

The above procedure is straightforward to implement in our case. Dropping the dimen¬ 
sionful prefactor in Eq. ([^, we write the differential distribution in terms of the Dalitz 
density as follows. 


pit) = I]ci/i(t), 

i 


(35) 
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where 


fi{t) = Po{t) +Pi{t), 
f2{t) = -7——rPo(^), 


{ruB 


fsit) = 


mlt^ 


{ruB - mo) 


mo) 

2Po{t), 


and 


Cl — 1, C2 — an, C3 — 


The dimensionless parameter an is dehned as 


an = 


"m^tan/dy 


rriH 


(36) 

(37) 

(38) 


(39) 


(40) 


We may now use the machinery of the optimized weighting procedure in order to calculate 
the statistical errors for a given model (i.e., for a given value of tan (3 /niH) and a given 
number of experimental events. We take as input some value of tan(3/mH, calculate the 
elements of the matrix M and then perform the sum in Eq. (^). The c° in this expression 
are given by the input values themselves and the Cj are as indicated^in Eq. ([39|). is fhns 
a polynomial in even powers of tan P/mH and is zero at the input value. The na error for 
a given experiment is simply gotten by setting x^=rP 


For our numerical analysis we have used the simple heavy quark symmetry forms for Eq, 
Fi and Fs (see Eqs. (p!^ ) - (|^) ), again taking ^{w) = 1 — 0.75 x (w — 1). In Fig. ^ we plot 
as a function of tan P/tuh for the case in which the input model is the SM {tan P /mH = 0). 
The solid, dashed and short-dashed curves correspond to iV = 10^, 3 x 10^ and 10^ events, 
respectively. The dashed vertical lines indicate the 2a upper bounds which could be placed 
on tan P/mH in each case. An approximate formula for this upper bound may be obtained 
from the expression for by setting y^=4 and truncating the polynomial in tan P/mn 
at the leading term. This leads to the following approximate expression for the 2a upper 
bound: 


n 30 

tanp/mn < GeVh 
For the three cases shown in the plot, this corresponds to 


(41) 


tan P/mH < 


0.039 GeV-\ 
0.052 GeV-\ 
0.069 GeV“\ 


N = 10^ 

AT = 3 X 10^ 
N = 10^ 


(42) 


in reasonable agreement with the exact results indicated in the plot. The upper bounds 
determined in this way would represent signihcant improvements on the current limits. 


®Note that we set ci = Cj* = 1, since we assume that the SM W-exchange contribution will always 
be present. 
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which are 0.46 GeV“^ and 0.38 GeV“^ (coming from the inclnsive semi-tanonic and tanonic 
B decays, respectively). Recall, fnrthermore, that the cnrrent theoretical nncertainties in 
the form factors would limit the reach of even an ideal experiment (with inhnite statistics) 
to about 0.06 GeV“^ (see Fig. ||). We see from Eq. (|4^) that already with about 3 x 10^ 
events one will have reached this limit. 

This approach is also extremely sensitive to non-zero input values of tanP/mH- As an 
example, we have calculated ^ 2HDM with tanP/mn = 0.1 GeV“^. The resulting 

plot, given as a function of tanP/mn, is shown in Fig. ||. The solid, dashed and short- 
dashed curves correspond again to = 10^, 3 x 10^ and 10^ events, respectively. The 
intercepts of the three curves at tan P/mu = 0 determine how well one can differentiate 
between this 2HDM and the SM in each case. Thus, for example, with about 10^ events one 
can differentiate between this model and the SM at approximately the 4.6(T level. 


V. CONCLUDING REMARKS 

In this paper we have examined the sensitivity of the exclusive decay B DtV to the 
tree-level charged Higgs contribution which is generic to type II two Higgs doublet models. 
The distribution in this decay is extremely sensitive to the ratio tan P/mu and may be 
used to appreciably improve the existing upper bounds on this quantity or, as the case 
may be, to measure a non-zero value. We have shown that while the existing theoretical 
uncertainties on this distribution are rather large, they will be reduced signihcantly once 
the analogous distribution for B —> D{e,^)v is measured more precisely. We estimate that, 
barring any further theoretical reductions in the uncertainty, this would eventually (i.e., 
assuming an inhnite number of experimental events) allow one to rule out a 2HDM with 
tanP/rriH > 0.06 GeV“^. We have also applied the optimized weighting technique to the 

distribution in order to calculate the minimum statistical uncertainty which could be 
attained for an experiment with a given number of events. The results of this analysis are 
very encouraging, showing, for example, that with just 10^ events one could rule out a 2HDM 
with tanP/rriH > 0.07 GeV“^. With the same number of events one could differentiate 
between a 2HDM with tanP/mu = 0.1 GeV“^ and the SM at the 4.6(j level. 

In the present work we have not made any use of the spin of the hnal state r. In principle, 
this extra observable could be used to improve the sensitivity of a given experiment. In 
practice, one has to allow the r to decay and study the distribution of its decay products. 
We have, in fact, used the optimized weighting procedure to study these distributions in the 
hadronic decay channels r — > ttz/ and r —> pv. The differential distributions in these cases 
may be written as functions of g^. Eh and cos 9h, where Eh and 9h are the energy and angle 
(with respect to the momentum of the D) of the meson h in the t—V rest-frame. (Recall 
that we wish to avoid observables which depend explicitly on the momentum of the r.) The 
upshot of this analysis is that for a hxed number of events, the upper limit on tan P/mn is 
improved by at most about 5%. In general, however, one can also expect a reduction in the 
number of events since one is now considering a very specihc decay mode of the r. Thus one 
is likely not to gain anything at all. It is thus our opinion that the g^ distribution represents 
perhaps the best tool which may be used in order to search for a charged Higgs signal in 
B —> DtV. 
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FIG. 



. Quark-level diagrams for the transition b —> ctu in a two-Higgs doublet model. 
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FIG. 2. Spin-0 and spin-1 contributions to the differential distribution in the standard model 
and in two Higgs doublet models with various values of tan (5/niH- The shaded bands labeled 
a, b, c and d correspond to the spin-0 contribution, (1 -|- po(t), with tan jd/mn = 0 (the 

SM), 0.06 GeV“^, 0.25 GeV“^ and 0.35 GeV“^, respectively. The solid curve corresponds to 
the spin-1 contribution, pi{t). The shaded regions indicate the current theoretical uncertainty 
in (1 -|- po{t) (due to the uncertainty in ({t) = Fo(t)/Fi(t)) if pi{t) is known exactly. 

Uncertainties in Fs{t) have not been included in this plot. 
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FIG. 3. Plot of the integrated width (normalized to G'|^|14{,pm^/1287r^) as a function of 
tan (5/niH ■ As in Fig. the shaded region corresponds to the theoretical uncertainty in 
C(t) = Fo(t)/Fi(t) and does not take into account the uncertainty in the scalar form factor, Fs{t). 
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FIG. 4. for the optimized weighting technique as a function of tan/3/m//. The “input 
model” in this case is the standard model (i.e., tan/3/m//=0). The solid, dashed and short-dashed 
curves correspond to 10^, 3 x 10^ and 10^ events, respectively. The vertical dashed lines indicate 
the 2(7 upper bound on tan/3/m// in each case, with the regions to the right of the dashed lines 
being excluded. 
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FIG. 5. for the optimized weighting technique as a function of tan/3/m// in a two Higgs 
doublet model. The “input model” in this case has tan/3/m// = 0.1 GeV“^. The solid, dashed 
and short-dashed curves correspond to 10“^, 3 x 10^ and 10^ events, respectively. The intercept at 
tan/3/m// = 0 shows that for 10^ events one can distinguish this scenario from the standard model 
at approximately the 4.Oct level. 
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